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Abstract. This paper reviews two different uses of the continuous wavelet trans- 
form for modal identification purposes. The properties of the wavelet transform, 
mainly energetic, allow to emphasize or filter the main information within mea- 
sured signals and thus facilitate the modal parameter identification especially when 
mechanical systems exhibit modal coupling and/or relatively strong damping. 



1.1 Introduction 

The concept of wavelets in its present theoretical form was first introduced in 
the 1970s by Jean Morlet and the team of the Marseille Theoretical Physics 
Center in France working under the supervision of Alex Grossmann. Gross- 
mann and Morlet [1] developed then the geometrical formalism of the continu- 
ous wavelet transform. Wavelets analysis has since become a popular tool for 
engineers. It is today used by electrical engineers engaged for processing and 
analyzing non-stationary signals and it has also found applications in audio 
and video compression. Less than ten years ago, some researchers proposed 
the use of wavelet analysis for modal identification purposes. However, use of 
wavelet transforms in real world mechanical engineering applications remains 
limited certainly due to widespread ignorance of their properties. 
This paper presents two ways of using the continuous wavelet transform 
(CWT) for modal parameter identification according to the type of measured 
responses : free decay or frequency response functions (FRF). 
In 1997, Staszewski [5] and Ruzzene [3] started to use wavelet analysis with 
the Morlet wavelet for the processing of free decay responses, but the identi- 
fication of the mode shapes was not performed. Wavelet transform allows to 
reach the time variation of instantaneous amplitude and phase of each com- 
ponent within the signal and makes the identification procedure of modal 
parameters much easier. Compared to other mother wavelets, Morlet wavelet 
provides better energy localizing and higher frequency resolution but has a 
disadvantage: frequency-coordinate window shifts along frequency axis with 
scaling while other wavelet types only expands the window. To remedy it, 
Lardies et al. [I] and Slavic et al. [5] preferred the Gabor wavelet. Argoul et 
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al. |6|7j chose the Cauchy wavelet and Le et al. [8] established a complete 
modal identification procedure with improvements for numerical implemen- 
tation, especially for a correct choice of the time-frequency localization. Edge 
effects are seen by all the authors and some attempts to reduce their negative 
influence on modal parameter identification are presented in [5], [8] and [9]. 
More recently, Staszewski [10] , Bellizzi et al. [H] and Argoul et al. [6j adapted 
their identification procedure in order to process free responses of nonlin- 
ear systems. In [6], the authors proposed four instantaneous indicators, the 
discrepancy of which from linear case facilitates the detection and charac- 
terization of the non- linear behaviour of structures. In [11) . an identification 
procedure of the coupled non- linear modes is proposed and tested on different 
types of non-linear elastic dynamic systems. 

Typically, CWT is applied to time or spatial signals. Processing frequency 
signals is unusual, but a first application was proposed in [TJ] and [13] ■ Ar- 
goul showed in [TJ] that the weighted integral transform (WIT) previously 
introduced by Jezequel et al. [14j for modal identification purposes can be 
expressed by means of a wavelet transform using a complex-valued mother 
wavelet close to the Cauchy wavelet. The WIT can be applied to either the 
ratio of the time derivative of the FRF over the FRF or the FRF itself. The 
representation with CWT provided a better understanding of the amplifica- 
tion effects of the WIT. When the FRF signal is strongly perturbed by noise, 
its derivative is hard to obtain; thus, Yin et al. [13] then proposed a slightly 
modified integral transform directly applied to the FRF and called the Sin- 
gularities Analysis Function (SAF) of order n. It allows the influence of the 
FRF's poles to be emphasized and direct estimation of the eigen frequencies, 
eigen modes, and modal damping ratios is performed from the study of the 
extrema of the WIT. 

1.2 Theoretical background for the continuous wavelet 
analysis 

A wavelet expansion uses translations and dilations of an analyzing function 
called the mother wavelet ip £ L 1 (M.) L 2 (M.) . For a continuous wavelet trans- 
form (CWT), the translation and dilation parameters: b and a respectively, 
vary continuously. In other words, the CWT uses shifted and scaled copies 
of ip(x) : ipb,a(x) — h^^^a^ whose L X (R) norms (H-H-J are independent of a. 
In the following, tp(x) is assumed to be a smooth function, whose modulus of 
Fourier transform is peaked at a particular frequency J?o called the " central" 
frequency. The variable x may represent either time or frequency; when nec- 
essary, the time and circular frequency variables will be referred to t and to 
ui respectively. The CWT of a function u{x) £ L 2 (M) can then be defined by 
the inner product between u(x) and ipb,a{x) 




(1.1) 
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where ip(.) is the complex conjugate of -0(0 • Using Parseval's identity, Eq. 
p.ip becomes 



T+[u](b i a) = — (u,$ b , a ) = — u(n)$(af2)e Mb df2 (1.2) 

where u{fi), tpb, a (^) and ip(f2) are respectively the Fourier transform (FT) 
of u(x), ipb,a(x) and ip(x); for instance for u(x) : u{Q) = u{x) e~ lQx dx. 
Moreover, when ip and u are continuous and piece- wise differentiable, and ip 
is square and absolutely integrable, and u is of finite energy, the CWT of u 
with ip is linked to the CWT of u with ip '■ 

T^,[u](b,a) = --T^[u](b,a) 

From p~5|, it can be seen that the CWT at point (b, fl = picks up 

information about u(x), mostly from the localization domain D(b, f2 = — ^) 
of the CWT, defined as 

D(b,Q = ^) = [b+ax, p ~aAx^,b+ax^+aAx, p }x[^~^ 1 ^ + ^±} (1.3) 
a a a a a 



where x^, and Ax^p are called centre and radius of ip, stated in terms of root 

1 r+oo 
Vila J — 00 



mean squares : x^p = .. hp J_°° x |"0(a;)| 2 dx and 



Ax^, = j^-y ( x ~ x i>) IV'WI dx, an d similar definitions hold on for 

the frequency centre Q^, and the radius AQ^ of ip. The area of D is constant 
and equal to four times the uncertainty : AAx^AQ^ = 4/x(0). The Heisenberg 
uncertainty principle states that this area has to be greater than 2. 
Referring to the conventional frequency analysis of constant-Q filters, the 
parameter Q defined as the ratio of the frequency centre Q$ to the frequency 
bandwidth {2AQ^); 

W 2 ao ± 2AH J [ ' 

a 

is introduced in [8] to compare different mother wavelets and to characterize 
the quality of the CWT. Q is independent of a. 

The notion of spectral density can be easily extended to the CWT (see [IB] 
for others properties such as linearity, admissibility and signal reconstruction, 
etc.). From the Parseval theorem applied to Eq. 11.21 it follows that 



+oo i r+ca 

ITvMM)! 2 db=— \u(k)\ 2 ip(ak) 



dk 



and it leads to the "energy conservation" property of the CWT expressing 
that if u£ i 2 (M), then Tj, e L 2 (R x R* + ,db^f), and ||T v ,[u]|| 2 = \\u\\ 2 = 
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u(JT) is hermitian iu(x) being real valued), one gets 



Tp^- ||w|| z . Finally, after changing the integration variable: a = ^ and because 



+CX) Z'+OO 



d6 



d/2 C, 

77 



7T ./,, 



+OG 



\u(f2)\ 2 df2 



r +oo |^(^)| 2 



where = J*^ duo < oo. 

This allows to define a local wavelet spectrum: E u< cwt{^, V) 



2 iT? 



(1.5) 



\T^[u]{b,%-) 



and a mean wavelet spectrum: E u _cwt{^) = E u cwt(^, b)db; and Eq. 



1.51 can be rewritten: J Q 



+ OC 



dfl = 0. 



When the mother wavelet ip is progressive (i.e. it belongs to the complex 
Hardy space i/)(f2) = : for : Q < 0), the CWT of a real- valued signal u is 
related to the CWT of its analytical signal Z u (see. [TB]) 



T^[u](b,a) = -Tji,[Z u ](b,a) 



(1.6) 



From [1731 the CWT has sharp frequency localization at low frequencies, and 
sharp time localization at high frequencies. Thanks to a set of mother wavelets 
depending on one(two) parameter(s), the desired time or frequency localiza- 
tion can be yet obtained by modifying its(their) value(s). In [8], two complex 
valued mother wavelets were analyzed : Gabor and Cauchy. Morlet wavelet 
is a complex sine wave localized with a Gaussian envelope and the Gabor 
wavelet function is a modified Morlet wavelet with a parameter controlling 
its shape. Cauchy wavelet i/j n of n order for n > 1, is defined by: 



n + l 



n+1 



Vx 2 + 1 



o % (n + l) Arctg{x) 



(1.7) 



The main characteristics of ip n (x) follow: its FT: ij) n (Q) 



2-irn" e 



■9(Q) 



where <9(i?) is the Heaviside function (ip n is progressive), its central fre- 



quency: £2o n = n, its L2 norm: HVVil^ = 



2 _ (2n)! 

*(n!) 



T 7r, its x- and frequency 



centres: = and 



its radius: A. 



I+jt , its frequency 



radius: AQ$ — " +1 , its uncertainty: 



1 + 2,1-1 > ^s quality factor: 



Q = ^2n+i and its admissi bility factor: C v , = 47r 2 ^ (2 ( ",^ )! = ^ H^nll^- 
The use of the Cauchy wavelet is legitimate when Q is less than 5 / as Q 
increases, both wavelets give close results, a little better with Morlet wavelet 
due to its excellent time-frequency localization, and behave similarly in the 
both time and frequency domains when Q tends toward infinity. 
Some numerical and practical aspects for the computation of the CWT are 
detailed in [8], especially when applied to signals with time bounded support. 
Let us note D u = [0, L] x [0, 2i^j^ yq \ the validity domain of the discretized sig- 
nal in the (6, (2) plane where $Nyq is the Nyquist frequency {fnyq 



2?^ with 
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sampling period T). The edge-effect problem due to the finite length (L) and 
to the discretization of measured data record and to the nature of the CWT 
(convolution product) is tackled. An extended domain D ext {b, fl) D D(b, fl) 
where f2 = ^-), is then proposed to take into account the decreasing prop- 
erties of tj) and ip by introducing two real positive coefficients c x and cq: 
D ext (b, fl = = [b+ ax^ - ac x Ax^,,b + ax^ + ac x Ax^ ] x \-± - cq—^- 

+ c fi—^-]- Forcing D ext to be included into D u leads the authors to de- 
fine a region in the (6, fl) plane where the edge effect can be neglected; this 
region is delimited by two hyperbolae whose equations are: fl = \c x Q^ and 
fl = j^c x Qij,^, and two horizontal lines whose equations are: fl — and 

fl = yq-^pr-- Let us consider a frequency flj for which u(fl) exhibits a peak 

and introduce a frequency discrepancy dflj (for example, the distance be- 
tween two successive peaks of u(fl)). From the intersection of the two hyper- 
bolae at the point (b = ^, fl = flj) and by imposing the frequency localiza- 
tion along the straight line: a = jf-, to be included into [flj — dflj, flj + dflj], 

some upper and lower bounds are found for Q: c^i^jj- < Q < -^P 1 and fi- 
nally, the authors proposed c x — cn — 5. 



1.3 Modal analysis and modal identification with CWT 

Experimental identification of structural dynamics models is usually based 
on the modal analysis approach. One basic assumption underlying modal 
analysis is that the behaviour of the structure is linear and time invariant 
during the test. Modal analysis and identification involve the theory of linear 
time-invariant conservative and non-conservative dynamical systems. In this 
theory, the normal modes are of fundamental importance because they allow 
to uncouple the governing equations of motion. Also, they can be used to 
evaluate the free or forced dynamic responses for arbitrary sets of initial 
conditions. Modal analysis of a structure is performed by making use of the 
principle of linear superposition that expresses the system response as a sum 
of modal responses. 

For linear MDOF systems with N degrees of freedom, the transfer function 
fflij (p) of receptance type is defined as the Laplace transform of the displace- 
ment at point j to the impulse unit applied at point i and it can be expressed 
as a sum of simple rational fractions 

»,<*>= f:(S+S) <-) 

where (^r)^ and are respectively the conjugate of the residues (A-)jj an d 
of the poles p r of (p). When the system is stable, p r = — £ r u r + iui r where 
u> r , uj r are the undamped and the damped vibration angular frequencies and 
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£ r is the damping ratio for the mode r. (Ar)- can also be expressed from the 



complex modes \r by = *' 1 where "f r = Xr ' GXr+^PrXr MXr, M. 

and C_ being respectively the mass and viscous damping matrices. Moreover, 
the FRF of receptance type Hij(uS) can be usually related to Mij(p) by: 
Hij(u)) — Mij(p — iuj), leading to 

XM-T.ii + , ) (") 

z — ^ V [OJ — UJ r ) — l£ r L0 r [u> + to r ) — l£ r UJ r J 

The free responses at point j in terms of displacements Uj (f) can be 
expressed according to the modal basis of complex modes: Xjr = Vjr + * K jr 
(r being the mode number, 1 < r < N) 



N 

u y~.) {t) = A rj (t) cos (* ri (t)) (1.10) 



where A r j(t) — \Xjr\ Pr e ^ rUJrt and $ r j(t) — ^rt ~ ^Pr + arctan ^^-J ■ Pr an d 

are defined from initial displacement and velocity of mode r (cf. |17|). 
In the case of proportional viscous damping (Basile assumption), introducing 
the real eigenvectors W r of mode r for the associated conservative system, and 

replacing the residues {A r )^ with ^ ~ in Eq. 11.91 lead to 
y ' m r (a>£ — cij 2 + 2i £ r w r w) 

r— 1 v ' ' 

Moreover, in Eq. 11.101 Hj r = 0; thus |xjr| can be replaced by (#r)j and the 

term arctan ^^-^ can be equal to or it according to the sign of the real 
part rjj r . 



1.3.1 Modal identification using free decay responses 

The processed signals Uj free) (t), whose expression is given in Eq. 11.101 can 
be considered under the assumption of weak damping (£ r <C -^j), as a sum 
of N modal components: A r j(t) cos(<5 r j (t)) consisting in asymptotic signals. 
The approximation of asymptotic signal means that the oscillations result- 
ing from the phase term: <P r j(t) are much faster than the variation coming 
from the amplitude term: A r j(t); it entails that the analytical signal as- 
sociated with A r j (t) e l * rj W can be approximated by: A r j (t) e l * rj W . There- 
fore, Eq. {OJ and the linearity of the CWT entail that: [u^ ree) ] (6, a) = 

5 S^Li Tip[A r j(t) e l<s>rj ^](b, a). Taking the Taylor expansion of each ampli- 
tude term of this sum, it follows that 

, N 

2V[uf" 5 ](M) = (A rj (b)e^^$(aU; r ) + R rj (b,a)^j . (1.12) 

r=l 
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[abs(CWT) - channel 2] 




Fig. 1.1. Modulus of FT and of the CWT for the 4DoF system (Q = 18) 

Assuming that each R r j(b, a) is small enough to be neglected, \T^[u < j r " ) ](b, a)\ 
is peaked in the (6, a) plane near straight lines of equation: a = $rjr\ = S 2 -- 
When each modal pseudo-pulsation tv r is far enough from the others, this 
straight line's equation easily provides an approximation of ui r - Eq. 11.121 im- 
plies that the CWT of asymptotic signals has a tendency to concentrate near 
a series of curves in the time-frequency plane called ridges, which are directly 
linked to the amplitude and phase of each component of the measured sig- 
nal. In [8], a complete modal identification procedure for natural frequencies, 
viscous damping ratios and mode shapes, is given in the case of viscous pro- 
portional damping and applied to a numerical case consisting in the free decay 
responses of a mass-spring-damper system with four degrees of freedom (4- 
DoF). The procedure is applied here to the free responses of a 4-DoF system 
whose natural frequencies are: /i = 0.0984 Hz, f 2 = 0.1871 Hz, / 3 = 0.2575 
Hz and / 4 = 0.3027 Hz and modal damping ratios: £i = 0.0124, £ 2 = 0.0235, 
£3 = 0.0324, £ 4 = 0.0380. Fig. O shows the FT and the CWT of the dis- 
placement of the second mass for which the four eigenfrequencies are clearly 
visible. Fig. [T73] presents for the third mode and for the four masses, the mod- 
ulus of the CWT and its logarithm that will allow to estimate the value of the 
damping ratio £3 (cf. [8]). The edge effect is delimited by two hyperbolae in 
the time frequency plane; Q being chosen to 20 (15.3 < Q < 43.1). Identified 
values are very close to the exact ones and identification errors are negligible 
inside the domain D ext . This procedure was also applied in [7j to a set of 
accelerometric responses of modern buildings submitted to non destructive 
shocks. From the CWT of measured responses, the ridge and the correspond- 
ing amplitude and phase terms for the two first modes were extracted. For the 
first instantaneous frequency uii, the processing revealed a slight increase of 
ui\ just after the shock. The origin of this non linear effect was attributed to 
the non-linear behaviour of the soil-structure interaction. Other preliminary 
results can be found in [6] where instantaneous indicators based on CWT 
are computed from the accelerometric responses of a non-linear beam to an 
impact force. A Duffing non-linearity effect was then identified thanks to the 
first and the super-harmonic components. 
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Fig. 1.2. Modulus of FT and of CWT around the 3rd mode for all channels. Loga- 
rithm of the modulus of the CWT (L = 250s, T = 0.175s, c t = c n = 4.7, Q = 20). 



1.3.2 Modal identification using FRF functions 

The first definition of the SAF of order n (n € N) uses the derivative of a 
transfer function M(p) of the linear mechanical system, where p = a + ib : 



H„(&, a ) = (-l)"2™^^{H(p)} 



(1.13) 



H„ (b, a) can also be expressed as the FT of the impulse response function 
h(t) filtered by: 2na^~t n e~ at O(t) where 0(t) is the Heaviside function : 

H„(&,a) = [ t n e~ at h(t) e" lb * dt (1.14) 

Finally, H„(6, a) is equal to the CWT of the FRF H(oj) with a mother wavelet 
being the conjugate of the Cauchy wavelet of order n multiplied by "l t : 



Hn(M) = -^rTf [H](b,a) = — / H 
a 2 a 2 jr 



(u)tf> n ( - — - ) du (1.15) 



Let us now consider a transfer function restricted to only one term F\ iUo (u>) 
of the sum given in (II. 9|) and perturbed by a Gaussian white noise n{ui) so 
that : H(u) = F\ tU , (u>) + n{uj) — + n(w). Let us introduce the 

coefficient p 2 defined as the ratio of the square absolute value of the CWT 
of the pure signal upon the common variance of the CWT of the noise itself: 

P 2 = ETT^TTT^TlTi = irrrr 2 T~~Ti2i ; it looks like a signal to noise ratio. 
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The use of the Cauchy-Schwarz inequality gives that : p 2 < \\F\ jUm \\ and the 
maximum of p 2 is reached for (b, a) — (ljq,X) when ipb ta (u}) — K F a ,b(uj) — 
K UjZtis-ia where K is a constant. In conclusion, ip(w) = K B^j —KB ip (u>); 
the analyzing functions if)(u>) which allow the minimization of the noise effect 
on the peaks of \T^,[H](b, a)\ are proportional to ^ ( s ) defined in Eq. 11.71 
The above definitions of the SAF of order n (n £ N* ) allow to better under- 
stand its effect when applied to discrete causal linear mechanical systems. In 
Eq. 11.131 the successive derivatives of a rational fraction make the degree of 
the denominator increase and cause an amplification of the pole's effect. So, 
when Eq. 11.81 is introduced in Eq. 11.131 the SAF of order n becomes 

U nM = 2,n\a^Y.{ A + f ) (1.16) 

r=l \ \ a + lb ~ Pr) (a + lb-Pr) ) 

The absolute value of H n (b, a) exhibits 2N maxima located symmetrically 
from the vertical axis in the phase plane. The effect of parameter n is to facil- 
itate the modal identification, especially when the structure has neighboring 
poles. As soon as n is large enough, the coordinates (b max , a max ) of these 27V 
maxima allow the estimation of the real and imaginary parts, respectively 
5i {p r } and 9 {p r }, of each of the 2N conjugate poles of the system : 

K {p r } = -£r w r w -a maXr and 3 {p r } = uj r w 6 maXr (1.17) 

The residue A r can be then estimated from the complex valued amplitude of 
the extremum called H™ aXr = H™ ax (a maXr ., 6 maXr ) 

n+l 

Ar w 2 "( fl max J 2 Xr (1 lg) 

7rn! 

As n increases, the estimation of the extrema is better for signals without 
noise; but for noisy signals, noise effects are also increased. Recently, Yin et al. 
[18] proposed the use of analyzing wavelets which are the conjugate of Cauchy 
wavelets in which n is replaced by a positive real number. Applications both 
to numerical and experimental data showed the efficiency of this technique 
using SAF [Hj. It was applied to the FRFs obtained with HI estimators, of a 
test structure designed and build by ONERA for the GARTEUR-SM-AG19 
Group. This structure was made of two aluminum sub-structures simulating 
wings/drum and fuselage/tail. Finally, the results obtained with SAF were 
very similar to those obtained with the broadband MIMO modal identifica- 
tion techniques provided by the MATLAB toolbox IDRC. 

In conclusion, the pre-processing of measured signals (free decay responses 
or FRFs) with CWT proved to be an efficient tool for modal identification. 
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